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ABSTRACT 



Studies of cool white dwarfs in the solar neighbourhood have placed a limit on 
the age of the Galactic disk of 8-9 billion years. However, determining their cooling 
1 ages requires the knowledge of their effective temperatures, masses, radii, and at- 

mospheric composition. So far, these parameters could only be inferred for a small 
number of ultracool white dwarfs for which an accurate distance is known, by fit- 
ting their spectral energy distributions (SEDs) in conjunction with a theoretical 
mass-radius relation. However, the mass-radius relation remains largely untested, and 
the derived cooling ages are hence mo del- dependent. Here we report direct mea- 
surements of the mass and radius of an ultracool white dwarf in the double-lined 
eclipsing binary SDSS J013851. 54-001621.6. We find M WD = 0.529 ± 0.010 M and 
-Rwd = 0.0131 ± 0.0003 R,0. Our measurements are consistent with the mass-radius 
relation and we determine a robust cooling age of 9.5 billion years for the 3570 K 
white dwarf. We find that the mass and radius of the low mass companion star, 
M sec = 0.132 ± 0.003 M Q and R scc = 0.165 ± 0.001 R Q , are in agreement with evolu- 
tionary models. We also find evidence that this > 9.5 Gyr old M5 star is still active, far 
beyond the activity lifetime for a star of its spectral type. This is likely caused by the 
high tidally-enforced rotation rate of the star. The companion star is close to filling its 
Roche lobe and the system will evolve into a cataclysmic variable in only 70 Myr. Our 
direct measurements demonstrate that this system can be used to calibrate ultracool 
white dwarf atmospheric models. 

Key words: binaries: eclipsing - stars: fundamental parameters - stars: late-type - 
white dwarfs 



1 INTRODUCTION 

White dwarfs are born hot and cool gradually over bil- 
lions of years. Their cooling is understood well enough to 
make them useful in meas uring the ages of stellar popula- 
tions (|Fontaine et aljfeOQlf) . White dwarfs with brown dwarf 
companions can be used to place constraints on the age 
of th e brown dwarf ijPinfield et alJ 120061 : iDav -Jones et al.l 
|200S| ). The growing number of wide field surveys, such as the 
UKIRT Infrared Deep Sky Survey (UKIDSS) and the Sloan 
Digital Sky Survey (SPSS) , have lead to an increase in the 
number of these systems (|Girven et alJ 1201 ll : ISteele et alJ 
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2011). A proper understanding of white dwarf cooling is es- 
sential for placing reliable limits on the ages of these systems. 

White dwarf temperatures, and hence ages, are deter- 
mined by fitting their spectral energy distributions (SEDs) 
using models of their atmospheres. At low temperatures 
(T c s < 6000 K) the atmospheric models need to in- 
clude the effects of collisions between hydrogen molecules 
l|Bergeron et al.lll995l ; fBergeron fc Leggettl 120021 ) (and with 
helium, if present). This effect dominates at near-infrared 
wavelengths in ultracool white dwarfs (T c g < 4000 K), sup- 
pressing the i nfrared flux and causi ng it to emerge at shorter 
wavelengths l|f3ergeron et al.| [T994l. In addition, the SEDs 
of ultracool white dwarfs depend on the white dwarf's sur- 
face gravity, and hence its mass and radius. Current white 
dwarf atmosphere models are not yet able to produce sat- 
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Table 1. Journal of photometric and spectroscopic observations. Exposure times for the X-shooter observations are for the UVB, VIS 
and NIR arms respectively. 



Date at 


Instrument 


Filtcr(s) 


Start 


Orbital 


Exposure 


Conditions 


start of run 






CUT) 


phase 


time (s) 


(Transparency, seeing) 


2011/11/01 


ULTRACAM 


ugr 


23:51 


0.90-2.25 


4.0 


Variable, ~ 1 arcsec 


2011/11/30 


RATCam 


r 


20:57 


0.81-1.07 


10.0 


Good, ~ 1.5 arcsec 


2011/11/30 


RATCam 


i 


22:41 


0.77-1.03 


10.0 


Good, ~ 1.5 arcsec 


2011/12/01 


RATCam 


r 


21:24 


0.80-1.08 


10.0 


Good, ~ 2 arcsec 


2011/12/02 


RATCam 


r 


23:40 


0.82-1.11 


10.0 


Good, ~ 2 arcsec 


2011/12/25 


X-shooter 




00:58 


0.85-2.25 


606,294,100 


Variable, ~ 1 arcsec 


2012/01/08 


RATCam 


i 


21:10 


0.79-1.05 


10.0 


Good, ~ 2 arcsec 


2012/01/14 


RATCam 


i 


20:34 


0.95-1.21 


10.0 


Good, ~ 1.5 arcsec 


2012/01/18 


ULTRACAM 


ugi 


19:43 


0.29-1.47 


4.0 


Good, 1.5 — 3.0 arcsec 



isfying fits to the observe d SEDs of ultracool white dwarfs 
l|Giammichele et alJlioi^ - ). 

If an accurate distance is known, the absolute magni- 
tude and a mass-radius relation can be used in conjunction 
with the SED modeling to estimate the mass and the surface 
gravity of the star. However, parallaxes are available only for 
a small handful of ultracool white dwarfs. In one of the best- 
studied cases, LHS 3250, this method gives an unrealistically 
low surface gravity of logp = 7.27, underlining the u ncer- 
tainti es in the SED mo dels l|Bergeron fc Leggettll2002h . Re- 
centlv lKilic et all (|2012T ) used this approach to determine the 
cooling ages of the ultracool white dwarfs SDSS J1102+4113 
and WD 0346+246, finding cooling ages of 10+°^ Gyr and 
11.2 ^°'g Gyr respectively. However, they still had to rely on 
the mass-radius relationship in order to determine all of their 
parameters. The majority of the ultracool white dwarfs have 
no parallaxes, and a canonical surface gravity of logg = 8 
was assumed for their analysis. However, altering this by a 
plausible ±0.5 dex changes t he resulting cooling age by sev- 
eral Gyr l|Kilic et al.ll201o"bT ). Furthermore, the mass-radius 
relation for cool white dwarfs is all but untested observation- 
ally, further adding to the uncertainties. Hence, at present 
no ultracool white dwarf has a reliable mass determination, 
and hence their cooling ages are subject to large uncertain- 
ties. 

One exception to this is the study of white dwarfs in star 
clusters. In the best cases the total age of the stars is well 
known as is the mass at the turnoff, hence the mass at the 
tip of the white dw arf cooling sequence can be measured 
(Hans en et alj [20071. Howev er, this picture is complicated 
by the presence o f binaries dBedin et al.ll200St) and helium 
core white dwarfs (|Kalirai et al.ll2007l ), which add additional 
complexity to the white dwarf luminosity function. Further- 
more, these white dwarfs cannot be used to constrain the 
age of the Galactic disk since they originated from a differ- 
ent population of stars. 

Double-lined eclipsing binaries can be used to mea- 
sure masses and radii with very few a ssumption s and 
to accuracies of better than 1 per cent (|Andersenlll99ll : 
ISouthworth et all l2007l l . independent of model atmosphere 
calculations. This method has been applied to white 
dwarfs in eclipsing binaries and has resulted in the 
most precisely measured white dwarf masses and radii to 
date jParsons et alJlibTol ; IPvrzas et aill2012i ; [ Parsons et al.l 
l2012ah . The subject of this paper, SDSS J013851.54- 



001621.6 (henceforth SDSS 0138-0016) was one of a number 
of candidate eclipsing white dwarf plus main sequence bina- 
ries identified in the multi-epoch SDSS photometric survey, 
known as Stripe 82 l|Becker et al.ll201ll 1. 

Here we present high precision photometry and spec- 
troscopy of SDSS 0138-0016, confirming its binary and 
eclipsing nature, and use these data to directly measure the 
masses and radii of both stars in the binary. The same data 
also yield the temperature of the white dwarf and hence the 
age of the system. 



2 OBSERVATIONS AND THEIR REDUCTION 

2.1 WHT /ULTRACAM photometry 

SDSS 0138-0016 was observed with ULTRACAM mounted 
as a visitor instrument on the 4.2m William Herschel Tele- 
scope (WHT) on the 1st November 2011 and the 18th Jan- 
uary 201 2. ULTRACAM is a high-speed, triple-beam CCD 
camera l|Dhillon et all l2007l l which can acquire simultane- 
ous images in three different bands; for our observations we 
used the SDSS u, g and either r or i filters. A complete log 
of these observations is given in Table [TJ We windowed the 
CCD in order to achieve exposure times of ~ 4 seconds which 
we varied to account for the conditions; the dead time be- 
tween exposures was ~ 25 ms. It is also possible to increase 
the exposure time of the u band observations by coadding 
the exposures on the CCD before readout. Since SDSS 0138- 
0016 is faint in this band (it' = 20.55), we used 5 coadds for 
the u band observations, resulting in exposure times of ~ 20 
seconds. 

All of these data were reduced using the ULTRACAM 
pipeline software. Debiassing, natfielding and sky back- 
ground subtraction were performed in the standard way. 
The source flux was determined with aperture photome- 
try using a variable aperture, whereby the radius of the 
aperture is scaled according to the full width at half maxi- 
mum (FWHM). Variations in observing conditions were ac- 
counted for by determining the flux relative to several com- 
parison stars in the field of view. The data were flux cal- 
ibrated by determining atmospheric extinction coefficients 
in each of the bands in which we observed and we calcu- 
lated the absol ute flux of our tar get using observations of 
standard stars (|Smith et al.l I2002T ) taken in twilight. Using 
our extinction coefficients we extrapolated all fluxes to an 
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Figure 1. Averaged X-shooter spectrum of SDSS 0138-0016, no telluric correction has been applied. Grey regions cover areas of severe 
(light grey) or nearly complete (dark grey) telluric absorption. The red points are the SDSS ugriz and 2MASS JHK magnitudes. 



airmass of 0. The systematic error introduced by our flux 
calibration is < 0.1 mag in all bands. 



2.2 LT/RATCam photometry 

Six primary eclipses of SDSS 0138-0016 were obtained in the 
r and i bands (three in each band) using RATCam, an op- 
tical CCD camera mounted on the robotic 2m Liverpool 
Telescope (St eele et al.ll20o"3 ). Each eclipse observation was 
composed of 75 x 10 second exposures. We used 2x2 bin- 
ning resulting in a readout time of ~ 5 seconds between 
exposures. These observations are summarised in Table [T] 

The raw data are automatically run through a pipeline 
that debiases, removes a scaled dark frame and flat-fields the 
data. The source flux was determined with aperture photom- 
etry using the ULTRACAM pipeline. The same nearby stars 
used to flux calibrate the ULTRACAM data were used to 
calibrate the RATCam data. 



2.3 VLT/X-shooter spectroscopy 

We obtained serv ice mode observations of SDSS 0138-0016 
with X-shooter jD'Odorico et alj l2006h mounted at the 
VLT-UT2 telescope. The observations were designed to 
cover an entire orbit of the system. Details of these obser- 
vations are listed in Table [1] X-shooter is a medium reso- 
lution spectrograph consisting of 3 independent arms that 
give simultaneous spectra longward of the atmospheric cut- 
off (0.3 microns) in the UV (the "UVB" arm), optical (the 
"VIS" arm) and up to 2.5 microns in the near-infrared (the 
"NIR'arm). We used slit widths of 1.0", 0.9" and 0.9" in 
X-shooter's three arms and binned by a factor of two in the 
dispersion direction in the UVB and VIS arms resulting in 
a spectral resolution of 2500-3500 across the entire spectral 
range. 

The reduction of the raw frames was conducted us- 
ing the standard pipeline release of the X-shooter Common 



Table 2. Identified white dwarf emission lines, twd is the sys- 
temic velocity of the line and K^j is the measured radial velocity 
of the line. 



Line 


Wavelength 


Twd 


^WD 




(A) 


(kms' 1 ) 


(kms -1 ) 


H14 


3721.948 


108.5 ±3.9 


82.1 ±4.2 


H13 


3734.372 


115.6 ±7.3 


74.6 ± 9.5 


Fel 


3737.131 


106.5 ±3.3 


87.0 ±4.6 


Fel 


3745.899 


108.6 ±2.8 


85.3 ±4.1 


H12 


3750.152 


109.1 ±5.3 


83.6 ± 7.7 


Hll 


3770.634 


104.4 ±3.8 


84.4 ±5.2 


H10 


3797.910 


105.7 ±6.0 


90.4 ±8.8 


H9 


3835.397 


111.7 ±7.6 


78.8 ± 9.0 


H8 


3889.055 


100.4 ±1.5 


81.5 ± 2.5 


Can 


3933.663 


105.0 ±1.0 


86.5 ± 1.0 


Can 


3968.469 


106.6 ±1.8 


85.0 ±2.8 


He 


3970.074 


103.8 ±1.1 


87.0 ± 1.3 


m 


4101.735 


106.8 ±1.1 


87.5 ± 1.3 


Cal 


4226.728 


101.8 ±8.3 


80.4 ±4.9 


H-y 


4340.465 


106.4 ±1.0 


87.3 ± 1.0 


Fel 


4383.545 


104.6 ±1.8 


84.0 ±2.6 


H/3 


4861.327 


106.5 ±0.5 


87.0 ±0.6 


Mgi 


5167.322 


110.6 ±6.6 


88.8 ± 7.9 


Mgi 


5172.684 


105.1 ± 1.7 


87.1 ± 2.6 


Mgi 


5183.604 


102.0 ±1.2 


85.9 ± 1.8 


Ha 


6562.760 


104.5 ±0.3 


85.5 ± 0.4 



Pipeline Library (CPL) recipes (version 1.3.7) within ES- 
ORex, the ESO Recipe Execution Tool, version 3.9.0. The 
standard recipes were used to optimally extract and wave- 
length calibrate each spectrum. The instrumental response 
was removed by observing the spectrophotometric standard 
star Feige 110 and dividing it by a flux table of the same star 
to produce the response function. The wavelength scale was 
also heliocentrically corrected. 
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Figure 2. A trailed spectrum of the Call 3934 A line. Emission 
can be seen from both components. The strongest component 
is from the white dwarf's chromosphere whilst the weaker (but 
larger amplitude) component originates from the main-sequence 
star. The red lines (online version only) track the motion of each 
component. Due to the short duration of the eclipse, the spectra 
taken at phase zero still covered the out-of-eclipse phase hence we 
are unable to say whether the white dwarf emission component 
is eclipsed. There was a 10 minute gap in observations around 
phase 0.5, which was made for calibration reasons. 



3 RESULTS 

3.1 Radial Velocities 

Figure [T] shows the average spectrum for SDSS 0138-0016. 
The M star features dominate the spectrum, but there are 
also several emission lines that move in anti-phase to the 
absorption features of the M star. These emission lines orig- 
inate from the white dwarf's chromosphere as a result of 
accretion of material from the wind of the M star. They 
have been seen in other close white dwarf plus main se- 
quence Wiwies_andrehabjy track the motion of the white 
dwarf l|Tappert et al.|[2011 a b). A list of the unambiguously 
detected emission lines from the white dwarf is given in Ta- 
ble [2] though there are likely to be additional lines at longer 
wavelengths which are obscured by the dominant M star. A 
trailed spectrum of the Ca II 3934 A line is shown in Figure^ 
Each emission line was fitted with a combination of a 
first order polynomial and a Gaussian component. For all 
the Balmer lines and the Can lines there is also an emis- 
sion component from the M star due to activity, in all cases 
it is the emission component from the white dwarf that is 
stronger (see Figure [5] for example). For these lines we fit 
both components simultaneously using a combination of a 
polynomial and two Gaussians. For the M star, we fit the 
K I 7700 A absorption line. Although fits to other M star 
absorption features give consistent results, the Kl 7700 A 
line is the cleanest feature available. Figure [3] shows sinu- 
soidal fits to the measured radial velocities for both the 
white dwarf and the M star. Table [2] lists the fitted radial 
velocities for the white dwarf emission lines. We measure 
Ki = 86.5 ± l.Okms^and K 2 = 346.7 ± 2.3 km s -1 . The 
two radial velocity amplitudes give the mass ratio of the 
two stars q = M scc /M W d = 0.249 ± 0.003. 




0.04 0.06 0.08 0.1 0.12 0.14 



MJD(BTDB) - 55920 (days) 

Figure 3. Radial velocity fits to the Can 3934 A emission line 
from the white dwarf (open circles) and the K I 7700 A absorption 
line from the M star (filled circles). 



3.2 Light curve model fitting 

Figure [4] shows our light curves of SDSS 0138-0016 around 
the expected time of the eclipse of the white dwarf by the 
M star. Our data confirm the eclipsing nature of the sys- 
tem. The reduced depth of the eclipse at longer wavelengths 
confirms that the bluer white dwarf is being eclipsed. 

To measure the system parameters we fitted all the 
light curves using a code writte n for the general case of b i- 
naries containing white dwarfs (jCopperwheat et al.ll201(ih . 
It has been used i n the study of other white dwarf-main 
sequence binaries |Parsons et al.l |2010| ; iPvrzas et al.l |2012| ; 
IParsons et"aLll2012af T The program subdivides each star into 
small elements with a geometry fixed by its radius as mea- 
sured along the direction of centres towards the other star, 
Roche geometry distortion and beaming are also included. 
The code also calculates the white dwarf contribution to the 
overall flux. 

The parameters needed to define the model were: the 
mass ratio, q — M sec /MwD, the inclination, i, the stellar 
radii scaled by the orbital separation R scc /a and Rwd/o,, 
quadratic limb darkening coefficients for both the stars, the 
time of mid eclipse, To, the period, P and flux scaling factors 
for each star. 

The primary eclipse shape does not contain enough in- 
formation to determine the inclination and scaled radii of 
both stars simultaneously. However, the amplitude of the el- 
lipsoidal modulation is related to the mass ratio and Rs ec /a. 
Therefore, since we knew the mass ratio from our spectro- 
scopic observations we could use it as a prior constraint and 
hence we were able to measure these parameters simulta- 
neously. Unfortunately, this approach did not work for the 
u band data since the white dwarf dominates the overall 
flux in this band, suppressing the ellipsoidal modulation and 
making it much harder to fit. It is also of much lower signal- 
to-noise and may be affected by activity from the M star. 
Since the RATCam light curves only covered the primary 
eclipse they could not be used to determine accurate param- 
eters. However, we used them to measure the magnitudes of 
the two stars. We fitted all the RATCam light curves sep- 
arately using the parameters found from the ULTRACAM 
light curves and allowed the flux scaling factors to vary. 
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Figure 4. Primary eclipse light curves and fits. The data shown in the left hand panels were obtained using ULTRACAM, whilst the 
right hand panels are from RATCam. Two models are shown for each light curve, the best fit model and the same model but with no 
primary eclipse, demonstrating the relative contribution of each star in the different bands. 



For fitting the light curves, we phase folded the data 
and kept the period fixed as unity. The limb darkening of 
both stars was set using a 4-coefficient formula: 



4 

£ 

!=1 



(1) 



where [i is the cosine of the angle between the line of sight 
and the surface normal. For the secondary star, we use the 
coefficients for a T e ff = 2 900, logp = 5 main sequence star 
jClaret fc Bloemenll201ll ). For the white dwarf we calculated 
the 4 coefficients from a white dwarf model atmosphere with 
Twd = 3500 and log g — 7.9, folded through the appropriate 
filter profiles. We kept all limb darkening parameters fixed. 

We used the Markov Chain Monte Carlo (MCMC) 
method t o determine the distributions of our model pa- 
rameters (|Press et all 120071 '). The MCMC method involves 
making random jumps in the model parameters, with new 
models being accepted or rejected according to their proba- 
bility computed as a Bayesian posterior probability. In this 
instance this probability is driven by a combination of x 2 
and the prior probability from our mass ratio constraint. 
Table [3] lists the best fit parameters and their statistical 



Table 3. Light curve model parameters from Markov chain Monte 
Carlo minimisation. The limb darkening coefficients (ai), are also 
listed for each star. 



Parameter 


9 


r 


i 


i 


77.14° ± 0.04° 


77.19° ± 0.02° 


77.25° ± 0.07° 


r WD /a 


0.0200 ± 0.0007 


0.0205 ± 0.0005 


0.0212 ±0.0018 


''sec / & 


0.333 ±0.001 


0.330 ±0.001 


0.328 ± 0.001 


ai (WD) 


-0.129 


1.185 


1.499 


a 2 (WD) 


2.710 


-0.473 


-0.402 


a 3 (WD) 


-3.079 


-0.041 


-0.457 


a 4 (WD) 


1.161 


0.096 


0.305 


ai (sec) 


0.2083 


0.5288 


0.6659 


a2 (sec) 


1.1341 


0.2451 


0.5135 


«3 (sec) 


-0.4029 


0.4560 


-0.2498 


a4 (sec) 


0.0467 


-0.2531 


0.0307 



errors, along with the limb darkening coefficients used for 
both stars. The fits to the g, r and i band light curves all 
give consistent results. Figure [4] shows the fits to each band 
around the primary eclipse. We also show the best fitting 
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Orbital phase 

Figure 5. ULTRACAM r band full orbit light curve of 
SDSS 0138-0016. The red line (online version only) shows the best 
fit to the data. The out-of-eclipse variations are caused by the 
tidally distorted main-sequence star presenting a different sur- 
face area during the orbit. The amplitude of this modulation is 
related to the mass ratio (q = M scc /Mwd) an d the radius of the 
main-sequence star. Since we know the mass ratio from the spec- 
troscopy, the amplitude of the modulation, combined with the 
eclipse shape, allows us to measure the orbital inclination and 
radii of both stars. The lower panel shows the residuals to the fit, 
a small region of data was affected by clouds (around phase 1.6). 



models but with the primary eclipse turned off, to illustrate 
the contribution of each star in the various bands. 

The best-fit model to the full light curve is shown 
in Figure [SI We find an inclination of 77.19° ± 0.02° and 
a white dwarf mass and radius of 0.529 ± 0.010 Mq and 
0.0131 ± 0.0003 R© respectively. The surface gravity of the 
white dwarf is then log g = 7.926 ± 0.022. 



4 WHITE DWARF TEMPERATURE AND AGE 

Figure [6] shows the SDSS spectrum of SDSS 0138-0016 in 
black and that of a second eclipsing white dw arf plus main- 
seque nce binary SDSS 1210+3347 in gray (|Pvrzas et al.l 
|2012| ). These two systems are very similar in that they con- 
tain M5 main-sequence stars and cool white dwarfs. Despite 
the low temperature of 6000 ± 200 K for the white dwarf in 
SDSS 1210+3347 a substantial blue excess is still produced 
ijPvrzas et al.l[201^ 1 . No such excess is seen in the spectrum 
of SDSS 0138-0016 implying that the white dwarf is much 
cooler than the one in SDSS 1210+3347. 

From our light curve fits we measure white dwarf mag- 
nitudes of g = 20.242 ± 0.007, r = 19.079 ± 0.006 and 
i = 18.773 ± 0.020. We also measure u = 21.42 ± 0.06, how- 
ever, the u band magnitude is likely to be unreliable since 
we know that the M star is active and any activity on the M 
star will heavily affect the u band due to Balmer continuum 
emission, therefore we do not use the u band magnitude to 
constrain the temperature of the white dwarf. 

The colours of the white dwarf in SDSS 0138-0016 are 
shown in Figure [7] along with those of other cool and ultra- 
cool white dwarfs llHarris et al.lll999l, l200ll ; lHall et alJliooi 
iKilic et all bOlOd lbT lLeggett et all 120111 ) . We computed a 




tOOO 5000 6000 7000 8000 

Wavelength [A] 

Figure 6. SDSS spectra of SDSS 0138-0016 (black line) and 
the similar eclipsing white dwarf plus main-sequence binary 
SDSS 1210+3347 (gray line) which contains a 6000 K white dwarf. 
Both of these systems contain an M5 main sequence star and 
their spectra have been scaled to match in the i band. Therefore, 
any discrepancy at shorter wavelengths reflects differences in the 
white dwarf components. 



set of model atmospheres spanning a wide range in effec- 
tive temperatures and at mospheric He/H abundance ratios 
ijGiammichele et alJlibrl ). we also include the opacity from 
the red wing of La. In all cases the mass was kept fixed 
as 0.529 M Q . For low temperatures (< 4000 K), the evolu- 
tion of the cooling tracks in the r — i colour depends al- 
most exclusively on the temperature, and in g — r on the 
He/H abundance ratio. Ultracool white dwarfs become bluer 
in r — i with decreasing temperature, because of collision- 
ally induced absorption, and bluer in g — r with increasing 
He/H ratio. The colours of the white dwarf in SDSS 0138- 
0016 therefore unambiguously constrain both its tempera- 
ture, T = 3570 "^gQ K, and its atmospheric composition, 
log(He/H) = 0.3. The He abundance seems plausible, as 
we know that the white dwarf is accreting from the wind 
of its H-rich companion, and that the accreted material will 
be mixed within the de ep convection zone that is typical for 
very cool white dwarfs l|Dufour et al.|[2007l ) . The model that 
best fits the measured colours of the white dwarf gives a 
surface gravity consistent with the value obtained from our 
light curve fits. This internal consistency strongly supports 
the validity of our results and also independently confirms 
the accuracy of the mass-radius relationship of cool white 
dwarfs. Our final best-fit model to the temperature and he- 
lium abundance implies a cooling age of 9.5 Gyr for the 
white dwarf in SDSS 0138-0016, making it one of the oldest 
white dwarfs with an accurate cooling age. 



5 DISTANCE AND KINEMATICS 

We followed the prescription of iBeuermannl (|2006l ) to esti- 
mate the distance of SDSS 0138-0016. For M-dwarfs, the sur- 
face brightness near 7500 A, and depth of the Ti O band near 
7165 A are a strong function of the spectral typ e . IBeuermannl 
(2006) provides a calibration of the surface brightness Ftio 
defined as the the difference between the mean surface fluxes 
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Table 4. System parameter of SDSS 0138-0016. 



Figure 7. Colour-colour plot of cool white dwarfs. The dotted 
lines are cooling tracks for different atmospheric compositions, 
from top to bottom log(He/H) = 2.0,1.0,0.5,0.3,0.0,-1.0. In 
the calculation of these tracks, the mass of the white dwarf was 
fixed to 0.529 M , as determined from our spectroscopic and 
photometric fits. The best model atmosphere fit is found for 
log(He/H) = 0.3 and T = 3570 + _}™ K. Cooling ages and tem- 
peratures are given for the log(He/H) = 0.3, and we determine 
the age of the white dwarf in SDSS 0138-0106 to be 9.5 Gyr. 
Furthermore, the surface gravity for this solution is in agreement 
with our measured value, meaning that the radius of the white 
dwarf in SDSS 0138-0016 is completely consistent with current 
mass-radius relationships. 



in the bands 7450-7550 A and 7140-7190 A. Measuring the 
observed flux /no from the spectrum, the distance is then 
calculated as 



d- 



FtiO 
/tio 



(2) 



Given the high accuracy of i? aec determined from the 
light curve model, the main uncertainties in the distance 
estimate are the flux calibration of the spectroscopy, and 
the spectral type of the companion. Adopting a conserva- 
tive uncertainty in the spectral type of M5 ± 1.0, and us- 
ing the SDSS and the VLT/X-shooter spectrum, we find a 
distance of 52+" pc, with the uncertainty in the spectral 
type dominating the error balance. Since the shape of the 
M-dwarf in SDSS 0138-0016 is distorted by the white dwarf 
(see Section 16. ip we calculate the distance using the mea- 
sured radius in several directions (and the spectra taken at 
the corresponding phase at which that radius is visible). We 
find that the oblateness of the M dwarf has a minor affect on 
the measured distance, much smaller than the uncertainty 
in the spectral type. The absolute magnitude of the white 



dwarf implied by this distance is M g = 16.66 



+0.46 



which 



is in agreement with the best-fit SED model, M g — 16.89, 
further validating our results. 

SDSS 0138-0016 has a relatively large proper motion 
(p, a = 336.2 ± 3.8mas/yr, (is — 32.9 ± 3.8mas/yr) and 



Parameter 



Value 



Orbital period 

TO (MJD(BTDB)) 

Orbital separation 

Orbital inclination 

White dwarf mass 

White dwarf radius 

White dwarf log g 

White dwarf effective temperature 

White dwarf cooling age 

Secondary star mass 

Secondary star radius sub-stellar 

Secondary star radius polar 

Secondary star radius backside 

Secondary star radius side 

Secondary star radius volume-averaged 

Distance 



0.072 764 91(2) days 

55867.007405(6) 

0.639 ± 0.004 R 

77.19° ±0.02° 

0.529 ± 0.010 M 

0.0131 ±0.0003 R 

7.926 ±0.022 
hiio 



3570" 



-so 



K 



9.5 



+0.2 , 



-o.3 °y r 

0.132 ±0.003 M 
0.211 ± 0.001 R© 
0.157 ± 0.001 R 
0.183 ± 0.001 R 
0.163 ± 0.001 R 
0.165 ± 0.001 R 
52+" pc 



a systemic velocity, measured from the X-shooter spectra, 
of 84.5 ± 1.2kms _1 . Using these values and the measured 
distance gives space velocities of U = —99 ± 9kms _1 , 
V = +198±6kms _1 and W = -45±3kms~ 1 , which makes 
it likely that SPS S 0138-0016 is a member of the thick disk 
|Pauh et al.ll2006l ). 



6 DISCUSSION 

6.1 System parameters 

A full list of the physical parameters of SDSS 0138-0016 is 
given in Table [4] The secondary star's shape is highly dis- 
torted due to the presence of the nearby white dwarf. It fills 
91% of its Roche lobe (as measured towards the LI point). 
Therefore Table [4] lists the radius of the secondary star in 
various directions. For our final discussions we adopt the 
volume-averaged radii. 

Figure [8] show the mass-radius plot for the white dwarf 
in SDSS 0138-0016. Also plotted are a number of theo retical 
mass-radius tracks from iBenvenuto fc Althaus! fl999). Our 
measurements are in excellent agreement with the models, 
although not precise enough to distinguish between the dif- 
ferent hydrogen layer masses. Nevertheless, this consistency 
reinforces our temperature and age measurement. 

The mass-radius plot for low-mass stars is shown in Fig- 
ure [§] The measured mass and radius of the low- mass star 
in SDSS 0138-0016 are consistent with evolutionary mod- 
els. Also shown are a number of other precise mass-radius 
measurements from low-mass stars that are in eclipsing bi- 
naries with white dwarfs. The precision of these measure- 
ments demonstrates the potential of these systems for test- 
ing low-mass stellar models. The number o f these systems 
has increased rapidly in t he last few years (ISteinfadt et aT 



200S| ; IPvrzas et al.ll2009l, |2012|; iNebot Gomez- Moran et al 



2009l ; lDrake et al.ll2009l . l2010l ; lLaw et al.ll201ll) making them 
a valuable resource for testing the mass-radius relationship 
for low-mass stars. 

As previously noted, the X-shooter spectra of 
SDSS 0138-0016 show emission components in the Balmer 
and Call lines originating from the M star due to activity 
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0.55 
0.65 




Figure 8. Mass-radius plot for the white dwarf in SDSS 0138- 
0016. The gray lines are t heoretical mass-radius tracks from 
iBenvenuto fc Althausl dl999Tl where the first number is the tem- 
perature of the white dwarf and the second number is the expo- 
nent of the hydrogen layer fraction. 



(in addition to the components from the white dwarf) . This 
indicates that the M star is still active, despite its age. The 
M5 star has an age of least 9.5 Gyr, likely more than 10 Gyr 
when the main-sequence l ifetime of the whi te dwarf progen- 
itor is taken into account. IWest et all (|2008l ) list the activity 
lifetime of an M5 star as 7.0 ±0.5 Gyr, substantially shorter 
than the age of the M star in SDSS 0138-0016. Therefore, 
it is likely that the tidally-induced rapid rotation of the M 
star keeps it active and makes it appear younger. 



Figure 9. Mass-radius plot for low mass stars. The mass and ra- 
dius values for the M star in SDSS 0138-0016 are s hown as well as 
other s from eclipsing whi t e dwarf binaries in re d j Parsons et al.l 
l20ld: IPvrzas et ail boil Iparsons et all l2012allbh. Oth e r mea - 
surements are fr om iKnigge^tTlT l l201ll) ; ICarter et al.l l l201lf) ; 
Ofir et al 1 d2012h. The solid line is the 8 Gyr isochrone from 
Baraffe et al.l 1119981 ) whilst the dashed line is a 5 Gyr model from 



Morales et al.l 1120101 ) which includes the effects of magnetic activ- 



ity. Also shown is a zoom in on the region around SDSS 0138-0016. 



70Myr at which point it will have an orbital period of 1.66 
hours. Due to the long angular momentum loss time-scale, 
systems of this type, so close to mass transfer, are likely to 
be rare. 



6.2 Evolution of the system 

We reconstruct the past and future e volution of SDSS 0138- 
0016 usin g the tools des c ribed in Schreiber fc Gansickd 
(|2003h and lZorotovic et all (|2010l . l201lf >. 

Assuming that the only mechanism of angular momen- 
tum loss from the system is via gravitational radiation then 
SDSS 0138-0016 emerged from the common envelope 9.5 Gyr 
ago with an orbital period of 5.28 hours. 

Fixing the common envelope efficiency to a — 0.25 re- 
sults in a mass of the white dwarf progenitor of 1.83M0. The 
evolutionary time scale of the white dwarf progenitor in this 
case would have been 1.63 Gyr, giving a total age of the sys- 
tem of 11.13 Gyr. Allowing for values of a between and 1 
but insisting that the evolutionary timescale of the progeni- 
tor is less than 4 Gyr (i.e. the system must be younger than 
13.5 Gyr), leads to a range of progenitor masses between 1.39 
and 2.OOM . 

A total age of ~llGyr makes it more likely that 
SDSS 0138-0016 is a member of the thick disk (as the space 
motion implies). This makes it a fairly old member of this 
population but consistent with previous kinematic studies 
that found that thick disk stars have a mean age of 10±2 Gyr 
|Feltzing fc Bensbvl I2009T ) . However, there is still some un- 
certainty as to whether the Galaxy has a thick-thin disk 
bi-modality l|Bovv et alj|2012h . 

The continuing loss of orbital angular momentum will 
lead SDSS 0138-0016 to become a cataclysmic variable in 



7 CONCLUSIONS 

Using high-precision photometric and spectroscopic data we 
measure the mass and radius of the ultracool white dwarf 
and low-mass star in the eclipsing binary SDSS 0138-0016. 
We use this information and the colour of the white dwarf 
to determine its atmospheric composition, temperature and 
age. We find that the white dwarf has a temperature of 
3570 "^g 1 ^ K and has been cooling for 9.5 Gyr. We also 
find that the mass and radius measurements for both the 
ultracool white dwarf and the low-mass star are consistent 
with evolutionary models. This supports the use of theoret- 
ical white dwarf mass-radius relationships when attempting 
to determine the temperature of ultracool white dwarfs us- 
ing SED fitting and parallax measurements. 

We find that the activity lifetime of the main-sequence 
star has been greatly extended due to being forced to rapidly 
rotate. The system is very close to Roche lobe overflow and 
will become a cataclysmic variable in only 70 Myr. 

The opacity from collisionally induced absorption from 
hydrogen in the ultracool white dwarf atmosphere is 
strongest in the near-infrared, making this wavelength range 
particularly sensitive to the temperature and atmospheric 
composition. Therefore, future measurements of the near- 
infrared magnitudes for the white dwarf will improve the 
precision and accuracy of the temperature and composition 
of the white dwarf. 
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